

my.data <- read.table("./source/getdata/analysis_traces.txt", header=TRUE, sep=" ")
#my.bandwidth <- c(0,0,0,0,0)
my.bandwidth <- c(0.02,0.03,0.04,0.05,0.01)
my.nb_obs <- nrow(my.data)
my.nb_var <- ncol(my.data)
my.nb_points <- 100
my.nb_iterations <- 500  
my.epsilon <- 0.0001
my.max.range = 2

# mydata cannot contain any 0 !!!!!
# otherwise the estimation is screwed because of the P[X_k  in (b_k, b_k + h]] in equation 1.3. The problem is the "(" with b_k not included.
# try with 2 variables X_1 and X_2 and the following values:
# (10, 0) and (0, 50)
# For Q( (0,0) , 50 ) we get H(0,0) - H(0, 50) - H(50, 0) + H(50, 50) = 0 - 0.5 - 0.5 + 1 = 0 !

my.data[my.data == 0] <- my.epsilon
